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We consider an interacting, one-dimensional Bose gas confined in a split trap, obtained by an har¬ 
monic potential with a localized barrier at its center. We address its quantum-transport properties 
through the study of dipolar oscillations, which are induced by a sudden quench of the position of 
the center of the trap. We find that the dipole-mode frequency strongly depends on the interaction 
strength between the particles, yielding information on the classical screening of the barrier and 
on its renormalization due to quantum fluctuations. Furthermore, we predict a parity effect which 
becomes most prominent in the strongly correlated regime. 
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The study of elementary excitations is a fundamental 
aspect of many-body theories. For neutral quantum flu¬ 
ids, these excitations at low energy correspond to sound 
waves for homogeneous systems and to inhomogeneous 
collective modes with discrete frequencies for confined 
ones. The analysis of the latter in ultracold quantum 
gases has been the subject of intense experimental Bi¬ 
ll and theoretical [9HH] activity in the last decades. 
A variety of different excitation modes has been char¬ 
acterized, as the monopole (breathing), dipole (slosh¬ 
ing), quadrupole, and scissor modes, to cite the best 
known. An unprecedented precision has been reached 
in the measurement of their frequencies, becoming one of 
the most reliable tests for theoretical models and tools 
to investigate many-body phases and beyond-mean field 
effects 0 nsHni One of the most interesting aspects 
of these excitations is that their frequencies depend on 
the microscopic properties of the system, yielding infor¬ 
mation e.g. on the equation of state or on its superfluid 
properties. 

The analysis of collective modes allows in particular 
to investigate the interplay between strong interactions 
and confinement in low-dimensional geometries [EM]. 
In this work, we show that by using a specific confining 
geometry, i.e. a localized barrier at the center of a quasi- 
onedimensional harmonic trap, one can directly access 
the effect of quantum fluctuations, which play a major 
role in low dimensions. Effectively one-dimensional sys¬ 
tems have been realized in ultracold-atoms experiments, 
by employing optical lattices to create arrays of tubes 
or by creating a single atomic waveguide on an atom 
chip [22U25] . Strongly correlated phases are more acces¬ 
sible in one-dimensional gases [2d, where interparticle 
interactions may be tuned by confinement-induced res¬ 
onances, and because, counterintuitively, in one dimen¬ 
sion the interactions become dominant in the low-density 
regime where particle losses and three-body recombina¬ 
tion effects are reduced. The demonstration of the pe¬ 
culiar fermionized Tonks-Girardeau phase constitutes a 
beautiful example [27, 28,. 


We focus on the dipolar excitation mode of a one¬ 
dimensional (ID) ultracold Bose gas, i.e. on a periodic 
oscillation of the center of mass of the atomic cloud. In 
ultracold-gas experiments, this sloshing mode can be ex¬ 
cited by a displacement of the center of the confining po¬ 
tential. A localized barrier can be created by microscope- 
focused laser beams [29: or by a light-sheet repulsive 
potential [3oi on]. For a purely harmonic potential, as 
predicted by Kohn’s theorem [32, [33], the dipole mode 
has the same frequency as the harmonic trap for arbi¬ 
trary interactions. In the presence of the barrier, Kohn’s 
theorem does not apply. In this work we show that the 
dipole mode displays an interaction-dependent frequency 
shift which allows to estimate directly the effective bar¬ 
rier strength seen by the fluid. We also find a surprising 
parity effect in the oscillation frequency, which becomes 
important in the strongly correlated phase and can be un¬ 
derstood in terms of fermionic rather than bosonic trans¬ 
port processes. The dipolar oscillation of the cloud real¬ 
izes in fact a specific type of quantum transport across 
the barrier. Quantum transport phenomena are more 
and more explored with ultracold atoms [341136 . 

System and Main Results .— We consider a tight 
atomic waveguide, containing N bosons of mass m at 
zero temperature, confined to a ID geometry. The 
bosons interact with each other via a contact interaction 
v(x — x') = g 5(x — x'), and are subjected to the harmonic 
confinement Vh(x) = ^mou^x 2 along the waveguide. The 
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FIG. 1. (Color online). Sketch of the potential along the ID 
wire. At t — 0 the potential is displaced by Ax to induce the 
dipole oscillation of the center of mass of the Bose gas. 
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FIG. 2. (Color online). Time evolution of the center of mass 
position xcm (t) for various values of barrier strength A = 
Uo/Ehlh, with — flbJhj and for interaction strength g = 
0. Solid black line: real-time evolution obtained from the 
numerical solution of the quench problem, red dashed line: 
estimate using the energy gaps of the many-body spectrum 
- a single frequency (E\ — Eo)/h in the first panel, and an 
additional frequency for the second and third panel. 


waveguide also contains a localized potential barrier at 
the center of the harmonic confinement, Vb(x) = UoS(x), 
which gives rise to a split trap for the particles. The 
corresponding Hamiltonian reads 

N h 2 *2 i N 

n = E ~2^id^ +Uo5 ^ + 2 m ^ x2j+ 2 S fi( x i~ x j) • 

3 = 1 3 3 , 1 =1 

(1) 

As sketched in Fig. |TJ starting from the ground state in 
the split trap, we study the dynamics following a sudden 
quench of the harmonic and barrier potential positions, 
both being displaced at time t = 0 + by a small amount 
Ax < with = yjh/muoh the characteristic har¬ 
monic oscillator length. This induces a collective oscil¬ 
lation of the gas of dipolar nature, and in particular the 
center of mass (CM) of the gas starts to oscillate period¬ 
ically in time with a main frequency around the new 
center of the trap. For specific interaction regimes (non¬ 
interacting, weakly, and infinitely strongly interacting) 
we have found the full time-resolved dynamical evolu¬ 
tion of the cloud following the quench, an example being 
illustrated in Fig. [2] 

In the absence of the barrier, when only the harmonic 
confinement is present, the oscillation is purely sinu¬ 
soidal: only the sloshing mode is excited by the quench, 
at the frequency of the harmonic confinement. This 
Kohn’s theorem result [32l 133] . specific to the harmonic 
confinement, is valid for arbitrary interaction between 
the particles, which can be understood as a change of 


reference frame into an accelerated one m- 

In the presence of the barrier (see again Fig. |2j, we 
observe the appearance of an additional harmonic com¬ 
ponent and a frequency shift of the dipolar oscillation 
with respect to the harmonic confining one. This fre¬ 
quency shift, which is the main topic of our analysis, 
increases with the barrier strength, and depends on the 
interaction regimes between the particles (see Fig. [3ja, 
b)). The frequency shift directly reflects the strength of 
the renormalized barrier seen by the fluid. The barrier is 
maximally reduced at intermediate interactions as a con¬ 
sequence of the competition of classical screening - occur¬ 
ring for weak interactions and increasing with the interac¬ 
tion strength, and barrier renormalization by the quan¬ 
tum fluctuations of the density - occurring for strong 
interactions and, being conjugate to phase fluctuations, 
decreasing for increasing interactions. Furthermore, we 
find that the frequency shift of the dipole mode in the 
strongly correlated regime depends on the particle num¬ 
ber being even or odd. This is particularly striking for a 
bosonic system, which does not display parity effects in 
other observables e.g. in the persistent currents [38, 35. 
When the barrier potential is not placed at the center of 
the harmonic trap, the parity effect is still present, but it 
is modulated by the position of the barrier, see Fig. [3|c) 
and m ■ Remarkably, we find also that signatures of this 
mesoscopic effect remains visible at finite temperatures, 
Fig. id). 

We detail below the steps and methods employed to 
obtain Figs. [2] and [3] 

Exact solutions .— In the two limiting cases of nonin¬ 
teracting (NI) and Tonks-Girardeau (TG) gas, we find 
an exact solution for the dynamical evolution of the 
gas. In the TG case, the solution is obtained by map¬ 
ping the system onto a gas of non-interacting fermions 
subjected to the same external (time-dependent) poten¬ 
tial 0D S3- In both cases, starting from the analyt¬ 
ical expression for the ground-state wavefunction 4>o <0 
before the quench, we perform a numerical evolution of 
the ground-state wave function in real time |\k[p 0 (£)) = 
exp(— iEL^t/h)^^ 0 ). The resulting time evolution of 
the position of the CM x CM (t) = f dx xn(x,t), where 
n(x,t) is the time-dependent gas density, is illustrated 
in Fig. [2] A Fourier analysis of x CM (t) shows that the 
center of mass motion is characterized by a main dipolar 
frequency, corresponding to the difference between the 
energy of the first excited many-body state and of the 
ground state of cjd = (-£q — Eo)/h, as well as by a 
higher harmonic with frequency (£3 — Eo)/h, in agree¬ 
ment with many-body perturbation theory [40] . 

In the following we shall focus on the dipolar fre¬ 
quency co>d and evaluate it in various interaction regimes. 
For the NI case, the TV-particle ground-state energy is 
given by Eq 1 = Neo, and the first excited one by 
E f 1 = (N — l)^o + £1, where Ej are the single particle 
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FIG. 3. (Color online), (a, b) Frequency of the dipole mode ua as function of the interaction strength g for Uo/Eh_£h = 1, 
in different regimes: (a) GP and LL solution for N — 50,51 (orange, cyan and purple lines respectively), the dotted lines at 
intermediate interactions are a guide to the eye in the regime beyond the validity of the GP and LL approaches; (b) ED solution 
for N = 4, 5 (blue and red respectively). The green dashed lines in the panels (a, b) correspond to the limiting cases of NI and 
TG gas (g —>• 0, oo). (c) Dependence of c^d on the position of the barrier d, for g —>• oo; calculated with the exact TG technique 
and corresponding LL for Ko = 1 (inset, with the LDA density, and gray lines, using the exact TG density), (d) Density plot 
of the Fourier transform of the center of mass oscillation |9fte[(Lc C M(k;)] | as a function of frequency and temperature, in the TG 
limit, for Uo/E^ih. = 1 and N = 4. 


energies, solutions of 


-^ d x + U 0 S{x) + 1 mulx 2 


1p n = n 


( 2 ) 


which can be expressed in terms of Hermite and Wittaker 
functions [40} 143 ] , In the TG case, the many-body en¬ 
ergy spectrum coincides with the one of a Fermi gas. In 
particular, the ground state energy is Eq G = £k 

and the first excited state one is Ej G = £ k + £n- 

Both NI and TG limits are indicated in Fig. [3[a, b). 

The TG solution allows to readily get physical in¬ 
sight on the parity effect. Using perturbation theory 
in the barrier strength we have huia = = 

H + (^l^bl^F) - (« , 0 G |'Hb|3'o G )> where n b = 

J2j=i UoS(xj). Using the explicit form of the ground- and 
first-excited many body wavefunction in the TG limit, 
'J TG (:ri, ...x N ) = ni< J - <fc <jvsign(a: J - - x k ) det[ip n (x e )\, 
where n = 0...N — 1 for the ground state and n = 
0...N — 2, N for the first excited state, we readily obtain 

hhj A = buj h + c 0 (|V’jv(o)| 2 - iV’jv-i(o)! 2 ). (3) 


The single-particle orbitals ^(0) vanish for k odd and 
are finite for k even m , therefore co>d is larger or smaller 
than cJh depending on the number of particles being even 
or odd. In analogy to transport phenomena in Fermi 
gases, for the strongly correlated (fermionized) bosons 
the dynamics is determined by the states at the (effec¬ 
tive) Fermi level. By using the analytical expression of 
'ipk we also find that in the TG regime the frequency shift 
|cc?d — | scales as 1/y/N for large N (see [40] for details 

as well as for the nonperturbative calculation in the TG 
regime). As we shall see below, this coincides with a spe¬ 
cial case of the LL prediction. If the barrier is placed at a 
distance d from the center of the harmonic trap, the TG 


solution allows to calculate exactly the dipole frequency 
by solving the modified Eq. (S2) with barrier potential 
Uo5(x — d ) [TOjE]. As shown in Fig. [3jc) the parity ef¬ 
fect is modulated, displaying an oscillating behavior as a 
function of d. 


Exact diagonalization (ED).— To cover all the interac¬ 
tion strength regimes, for small numbers of particles, we 
use a numerical method based on the exact diagonaliza¬ 
tion of the Hamiltonian 0 : We calculate the low-energy 
eigenspectrum of the many-body system, and obtain the 
dipole mode frequency as uj d = (Ei—Eq)/H. To represent 
the Hamiltonian ([!]) we have chosen the Wparticle basis 
built using the single-particle eigenfunctions of Eq. (S2), 
which are analytically known. The truncation of the 
Hilbert space to a number S of single-particle states is 
the only approximation performed with this technique. 
Since the dimension of the Hilbert space rapidly grows 
with S and the particle number N according to ( 5+ ^ _1 ), 
we are limited to small N |45|. The predictions of ED 
calculations are shown in Fig. [3jb). 


Mean field approach .— For larger values of the par¬ 
ticle number we adopt complementary approximate ap¬ 
proaches. In the regime of weak interactions, neglect¬ 
ing quantum fluctuations, we describe the fluid as a 
Bose-Einstein condensate [46] by means of the mean- 
field Gross-Pitaevskii (GP) equation: [— ^dx + Uod(x) + 
\rmjj\x 2 + g|4>| 2 ]<F = where <F(x) is the condensate 
wave function, normalized to the particle number N , and 
fi the chemical potential. We integrate this equation in 
imaginary time to find its ground state solution, and, 
subsequently, we evolve the ground state in real time with 
the shifted Hamiltonian, to calculate the time evolution 
of the position of the CM, and thus the frequency of the 
dipole mode in Fig. ia). 
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Luttinger liquid theory. — In the strongly interacting 
regime we take into account the effect of quantum fluctu¬ 
ations using the Luttinger liquid (LL) theory [47]. This 
is a a low-energy, quantum hydrodynamics description 
of the bosonic fluid, in terms of the canonically con¬ 
jugate fields 0 and 0 , corresponding to fluctuations of 
the density and phase, respectively. The presence of the 
smoothly varying harmonic potential is taken into ac¬ 
count within the local density approximation (LDA), i.e. 
expressing the equation of state of the inhomogeneous 
system as a functional of the local density n(x). The 
effective LL Hamiltonian for the inhomogeneous system 
reads [48, [49] 


o/L L 

n 0 


A 

2ir 



v s {x)K{x){d x (j){x)) 2 + ^^l(d x e(x)) 2 
K [x) 


(4) 

The local Luttinger parameter K(pc) and the sound 
velocity v s (x) depend on the microscopic interaction 
strength g , and are determined according to v s {x)K(x ) = 
fv7rn(x)/m , and v s (x)/K(x) = d n fi{n{x))/lm. We as¬ 
sume an equation of state of the form /i(n) = gn u [T9 . 49l: 
this connects the GP regime, where fi{n) = gn , and TG 
regime, where p(n) = (h 2 tt 2 / 2m)n 2 . The parameters g 
and v are obtained from the Bethe ansatz solution of the 
homogeneous Lieb-Liniger model @ 01 EDJ. For this model, 
the dependence on the interaction strength of the Lut¬ 
tinger parameters is known (see, e.g. m); let us recall 
in particular that in the TG limit Kq = 1 and v s = vp, 
the Fermi velocity of the fermionized Bose gas. 

The localized barrier potential induces a very rapid 
spatial variation of the confinement, hence cannot be 
taken into account in the LDA description. We treat 
its effect perturbatively in the limit of weak bar¬ 
rier strength. The barrier, located in full gener¬ 
ality at position d, yields a non-harmonic contribu¬ 
tion to the Luttinger liquid Hamiltonian of the form 
= J^° oQ dxU()S(x — d)p(x). Keeping only the low¬ 
est harmonics in the density field expansion p{pc) = 
[n(x) + d x e(x)/n] £t_oo e 2iW{x)+2iln dx'n(x'), we ob _ 
tain the dominant term ~ 2Uon(d) cos[20(d) + 

2 tt J_ oo dxn(x)\. Integrating out the high-energy modes 
of the field 0(x), we obtain a renormalization of the 
barrier strength by quantum fluctuations, such that the 
effective barrier is given by £/ eff = Uoe~ 2G ( d \ where 
G(x) is the local two point correlation function for 6(x). 
In the particular case where the barrier is at the cen¬ 
ter of the trap and for large particle numbers we find 
C/ eff = Uo 1 ^) , where k, = Kq^v/2, and a is a non- 
universal parameter dependent on the cut-off of the ef¬ 
fective LL theory m- Perturbation theory |40 finally 
yields 


where f(u) = (I + I)r 2 (4 + 4)2?+ 2 /r(^ + l). The inte- 
gral of the density gives rise to an oscillating term as 
a function of d shown in Fig. j3^c), which reduces to 
(—1)^ for d = 0. This explains the parity-dependent 
frequency shift of the dipole mode. We have used Eq.([5]) 
to estimate the interaction-dependent frequency shift of 
the dipole mode in Fig. |3ja). For repulsive interactions 
k ^ 1 , and the effective barrier strength is smaller than 
the bare one, decreasing as the interaction strength is 
decreased from infinite to intermediate values (k, -a- 1 for 
g oo, and n —» oc for g —> 0 ), which explains the 
approaching of uj d to at decreasing interactions. Fi¬ 
nally, we notice that in the TG limit one has k = 1 and 
n TG ( 0 ) = (1/4) v^F (in LDA), thereby recovering 
the |cjd — ^h| 1/y/N scaling behavior found with the 

TG exact solution. 

Experimental considerations. — The above analysis is 
readily extended to include both thermal effects and a 
finite width of the barrier. We have performed an ex¬ 
act finite-temperature calculation in the TG limit H3) 
see Fig. id)- Interestingly, the main features of the 
dipole-mode frequency shift remain visible at finite tem¬ 
peratures: for temperatures k^T > thermal fluctua¬ 
tions mix the characteristic zero-temperature frequencies 
of odd and even number of particles. Finally, we have 
also checked that a finite barrier width of the order of 
the interparticle spacing does not considerably affect the 
estimate of the frequency shift [52]. 

Conclusions. — As a prototype of quantum transport 
in an interacting ID Bose gas, we have studied the dipolar 
oscillations of the gas across a localized barrier, induced 
by a sudden shift of the trap center. We have found that 
the main frequency of the oscillation allows to determine 
the effective barrier strength seen by the fluid. The full 
quantum solution also displays a peculiar parity effect, 
due to the combination of fermionic transport properties 
for the correlated Bose gas and the harmonic trap geom¬ 
etry. Observation of this parity effect and of the shift of 
the dipole frequency with the interaction strength would 
provide a non-ambiguous evidence of the effect of quan¬ 
tum fluctuations. 
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SUPPLEMENTAL MATERIAL FOR: 

“Dipole mode of a strongly correlated one-dimensional Bose gas in a split trap: parity 

effect and barrier renormalization” 


I. Dipole mode frequency from many-body perturbation theory 


We consider an excitation of the center of mass (CM) of the system via a sudden displacement of the external 
harmonic and barrier potentials V ex t = YljLi \mu\x 2 - + UoS(xj) by an amount Ax at time t = 0. For t < 0 the 
system is assumed to be in its many-body ground state |^q < 0 ) °f Hamiltonian H t<0 before the quench. In accordance 
to perturbation theory, we decompose the initial ground state into the eigenstates of the system after the quench, 


l^ <0 > = |^ u > + Er=Ufcl*D- where |¥ t / u ) = E\ 
0 


>t ^0 i ^0 


) and c k = (v k 


t> o 


\n t<0 - n^°\% )/(K - 

E^°). Assuming the displacement small, Ax we have ~ + Ax ^2^=1 d Xj V^°(xj ) + 0(Ax 2 ), where 


d x V^°(xj) = muo^Xj + Uo5'(xj). The expectation value of the position of the center of mass as a function of time is 


given by x CM (t) = f dxxn(x, t), in terms of the density n(x, t) = f daq,..., dxjy ^2j —i 5(x — Xj) |4/q (aq, ...,xjvd) | 2 
Perturbation theory finally yields 


/ « N oo 

Axx / dxi,... ,da?jyy^(5(a:-a;j) y^c fe 2 5fte[Wg(a;i,... ,Xjv)^fe(a:i,.. .,x N ) 

J j =l fc=l 


^—i(E k —Eq (Si) 


From this expression we see that the periodic oscillation of the center of mass motion is decomposed in a Fourier 
series at frequencies (E^° — E^°)/h, where the weight of each component depends on the coefficient c k and the 
overlap integrals of the many-body wavefunctions 


In the absence of the barrier, i.e. Uo = 0, it follows from Kohn’s theorem that the only non vanishing matrix 
element c k is the km 1 one, as we have checked in our numerical solution. Therefore the center of mass evolves at a 
single frequency uj q = ( E p° — E^ 0 )/h, ie the dipole one, where in particular we have that uq = ^h- 


In the presence of the barrier, i.e. Uo > 0, the matrix elements c k for k > 1 are in general non vanishing and the 
dynamics of the center of mass will be determined by more frequencies (E^° — E^°)/h. However, for small Uo, the 
contribution of the higher energy states with k > 1, gets less and less important, the most important contribution 
being the dipolar one. The dependence of the matrix element c k on the quantum number k is shown in Fig. [4] below. 
The overlap integrals in Eq.(Sl) further decrease the weight of the higher-frequency components. 


In conclusion, even in the presence of the barrier, we obtain that the sudden quench of the center of the trap excites 
mainly the dipole mode. The second harmonics of the motion is also obtainable with this method and well agree with 
Fig. 2 of the main text. 




FIG. 4. Absolute value of the matrix element c k as a function of the quantum number k , for the same value of dimensionless 
barrier strength considered in the main text A = 1, and different interparticle interaction strengths g/E^ih = 0, 3, 30 (from left 
to right). 
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II. Details on the exact solution of the single-particle problem 
A. Exact solution for a centred barrier 


The one-body Schrodinger equation in the presence of harmonic and centred barrier potentials reads 

^Pn — ^n^Pn • 


Jt? 1 

~2~ d l + «*) + 2 mU * x2 


(S2) 


Scaling all the quantities in units of the harmonic oscillator level spacing and characteristic length 

= y/h/muj^ we obtain the Schrodinger equation in the reduced form 


-\d 2 x + \5{x) + I x 2 


VVi — t 


(S3) 


where A = Uo/E^i^. The odd solutions ip 2 n-hi of the harmonic oscillator without barrier (A = 0) are still solution of 
this equation with barrier, and having a node at the position of the barrier, they are independent of the strength of 
the barrier potential. The even solutions ip 2 n can be expressed in terms of Whittaker functions, after noticing that 


Eq. (S3) for x > 0 corresponds to the differential equation for the parabolic cylinder functions [431E2] • 

02 n(x) =AfD Sn (\x\) , 

, . (7 r 7 r \ (7 r tt \ 

D Sri (x) = cos (^- + -e n j Yi - sin + -e n J Y 2 , 

where Af is a normalization constant, and 

Fl _ r (i-hd e - 4 M q £ + i I X 2 

^2 4’ 2 ’ 


(S4) 


r(f 


Y 2 = - V4 t q/ e -V^,M + y £,ar ! 


^2 2 


3 3 


4 2 ’ 


Here, M(a, 6 , c) is the confluent hypergeometric function and T(a) is the complete T-function. 

Imposing the cusp condition at the position of the barrier, d x ip n ( 0 + ) — d x ip n (0~) = \7p n (0), we get the following 
expression for the eigenvalues e n \ 


r (l - \e n ) = _A 

r (1 - K) 2 


(S5) 


We observe that for A -A 0 the energy eigenvalues e n tend to the even eigenvalues of the harmonic oscillator without 
barrier. In the opposite limit A —)> oo they tend to the odd eigenvalues, that become doubly degenerate. 


B. Exact solution for a non-centred barrier 

The one-body Schrodinger equation in the presence of harmonic and non-centred barrier potentials reads 

h 2 1 

— -— d 2 + UoS(x — d) + -rmJ^x 2 
2m 2 

In this case the parity symmetry of the system is broken, and we can not distinguish between even and odd solutions. 
The solution can be written piecewise in the form: 

tpn(x) = AfiQ(d - x)D Sn (-x) + Af r Q(x - d)D Sn (x) , (S7) 

where 0(a;) is the Heaviside function, and D £n (x) are the parabolic cylinder functions defined in Eq. ( |S4| . The 
condition of continuity and the cusp condition at the position of the barrier lead to the trascendental equation for the 
energy eigenvalues m- 

(£n - 1/2) {D en _ t (-d)D en (d) + £> 6B _i(d)I? en (-d)) = V^A D £n (-d)D £n (d) . 


= Zn'tpn 


(S6) 


(S 8 ) 
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III. Details on the inhomogeneous Luttinger liquid description 


In order to diagonalize the inhomogeneous Luttinger liquid (LL) Hamiltonian in the presence of harmonic confine¬ 
ment, Eq. (4) of the main text, we have introduced, in analogy to what is done for the quantum harmonic oscillator, 
the bosonic conjugate fields bj and fet in the following way: 





x)e zu m + <p*(x)e~^ t b 


j) ; 


(S9) 


such that the canonical commutation relation [— 9(x)/tt, d x >4>(x')] = i5(x — x') and the orthogonality relation 
f d xipj(x)ipi(x) = Sjj imply that [bj, b)} = Sjj. Hamiltonian (4) of the main text is then diagonal in terms of 
the bj and bj fields, 


Ho L = V froj. 


j=o 


b ) h 3 


(S 10 ) 


provided that the modes’ wavefunctions ipj(x) satisfy the differential equation 

-U}jy/v s (x)K(x)ipj{x) = v s (x)K(x)d x d x (y/v a (x)K(x)ipj(x)) 

This differential equation has the form of a Sturm-Liouville problem, and the functions pj (x) can be chosen 
real. The spatial dependence of the parameters K and v s can be determined in LDA, assuming a dependence 
/i(n) = rjn lf , yielding K{pc) = K 0 (l — x 2 / R 2 ) 1 /"- 1 / 2 , and v s {pc) = v^y/l — x 2 /R 2 , where R is the Thomas- 
Fermi radius which is determined from the normalization condition N = f R R dxn(x), Vo = and 

Kq = yj{dii/dn)\TG/{d^/dn)\ x= {) = yjh 2 ir 2 n{d) 2 ~ 1 'jmr\v. In LDA, we get from the equation of state the follow¬ 
ing expression for the density: n(x) = [p(p — V(x))] 1 /", where V(x) = moo^x 2 /2 and (i = mu^R 2 / 2 . 

Eq. dS10|) has an analytical solution 


(Sll) 




{Uj/u h ) 2 = (j + 1)(1 + jv/2) 


where Cj(x ) is the Gegenbauer polynomial and T(a) is the complete T-function. We notice immediately that the 
lowest eigenvalue j = 0, that corresponds to the frequency ujq of the dipole mode (j = 0), is given by for any 
interaction strength, in agreement with Kohn’s theorem. 

Let us consider now how the barrier term of the Hamiltonian H LL = Hq L + ^b L ’ ^t> L = f^°oo dxUo5(x — d)p(x ), 
affects the frequency of the dipole mode of Hq L . Considering the LL density field expansion [48], and keeping only the 
lowest harmonics / = ±1 in p{pc) = [n(x) + d x Q(x)/ tt\ e 2tW ^ +2tl7r dx n(<x \ we obtain as the most dominant 

term T~0^ ^ 2Uon(d ) cos[2 6(d) + 27r d xn(x)\. Since we are interested only in the dipole excitation mode, if the 

barrier strength is small compared to the characteristic energy of the dipole mode (U eff /£h < frcjd), we can integrate 
out all the higher modes, j ^ 1. Considering the Fourier decomposition of the field 6(x,t) = e ~ luJjt 0j{x) : and 

taking the zero-temperature average over the vacuum of excitations, leads to a renormalization of the barrier strength 
U eff = CTq( 0| cos(2 0j(d)) |0). The LL being an effective low-energy field theory, the sum over the modes can not 

extend up to infinity, but should stop at a certain cut-off j c , that is intrinsic of the effective theory and that can not 
be determined within the theory itself. In our calculations we take it to be proportional to the Fermi energy [54] . 
thus j c is proportional to the number of particles j c = 7V/c, up to some numerical factor c of order 1 . The latter is 
fixed imposing the matching between the LL solution and the exact TG one in the infinitely strong interaction limit, 
Kq 1 . We define ( 0 | cos (2 Y^=\0j(d))\d) = exp[—2G(d)], where G(d) = (0|(X)^i 0j(d))) 2 \0). Using the mode 


m 


decompositions (S9), the correlation function takes the form G{d) = (ir/2)K{d)v s {d) (l/ Co ’i)l ( ^i(^)| 2 - Thus, 

general, the effective barrier strength for the lowest mode is given by 

U eff (d) = U 0 e~ 2G ^ . (S 12 ) 
















In the particular case d = 0, when the barrier is placed at the center of the harmonic confinement, and in the 
N 1 limit, the correlation function can be easily evaluated using the analytical solution of Eq. (S10), giving 
G(0) = /c(l/2) log(7V/a), where a is a numerical factor proportional to the cut-off parameter c. Using this result, we 
obtain the expression for the effective barrier strength for the dipole mode: 


= V„ (£)' 


(S13) 


where k = K(0)v s (0)/uJhR = Kq y 2 . 

Having integrated out the higher modes and consequently renormalized the barrier strength, we can rewrite the 
barrier term of the Hamiltonian in terms only of the field 6>o, as 2 U ef[ (d)n(d) cos[2#o(d) + 2tt dxn(x)]. 

Upon a Taylor expansion for small 0q to second order, substituting Eq. (S9) for the field and neglecting constant 
contributions to the Hamiltonian, we get the final, diagonal, expression: 


o/LL 

n b 


2n(d)U ef£ (d)7TV s (d)K(d)^pl(d)u: h 1 cos(27r J dx n{x)) (b^bo + H.c.^j 


This, after substituting the analytical expression of po given by Eq. (SI 1), leads us to Eq. (5) of the main text for 
the shift of the dipole mode frequency: 


= 


n(d)U^(d)K(d)v s (d)f(^l)T\^^2i^ 


huj^R 


r(f+i) 


cos 27r 


f 


dxn(x) 


(S14) 


which, for d = 0, simplifies to 


-Wh = (—l) N n(0)Uo K 


vi((bH)T\h+h)2 i+ * 


ra+1) 


(S15) 


IV. Finite-temperature center-of-mass oscillation spectrum 
A. Time evolution in the TG limit 

We obtain the oscillation spectrum at finite temperature (Fig 3(d) of the main text) by performing a Fourier analysis 
of the time evolution of the position of the CM x CM (t) = f dx xtit(x, t), where ut(x , t) = ^)| 2 - Here, 

f{e) is the Fermi distribution function at finite temperature for a state of energy £, Sj are the single particle energies 
determined in Section II, and Vb( x , 0 are th e time dependent single particle wavefunctions determined by numerically 
evolving in real time with the after-quench Hamiltonian the exact initial-state wavefunctions determined in Section II. 


B. Linear response theory 


The dynamics of the center of mass of the system can also be analyzed, in the limit of small oscillation amplitude, 
through linear response theory. This allows to obtain the time evolution of the particle density n(x, t) = (p(x, £)), where 
p is the density operator in second quantized form, in response to the perturbation operator H p = f dxV p (x,t)p(x) 
where V p (x,t) = 0(t)Axd x Vj^°(x), see Section I. Within the framework of linear response theory the corresponding 
evolution is then given by (p(x,t)) = (p 7 (x,t)) + f dx' J d t'x(x,x';t — t')V p (x', £'), where y(x, x'; £, £') = (1 /ih)Q(t — 
t')([p I (x, £), p T (x' : t')]), and p 7 (x, t) = e' int "° t ^ h p(x)e~' int "° t ^ h gives the unperturbed evolution of the density operator. 
Recalling that the position of the center of mass is given by x CM (t) = J dx x(p(x,t)) we readily obtain in Fourier 
space Sx cm (lj) = f dx x f dx' y(x, x'; u)V p (x', cj), with 


X(x,x’;uj) = ^ (rn\p(x)\n)(n\p(x')\m)e 


-I3E„ 


n^rri 


1 


(cj - (E n - E m )/h)) +i0+ 


(cu + (E n — E m )/h)) + i0+ 4 

(S16) 


Here, Z(/3) = e ^ m , (3 — l/k^T, n and m denote many-body states. 
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In the TG limit, using the Bose-Fermi mapping, the density-density response function coincides with the one of a 
non-interacting Fermi gas: 

*<„» = (i/B)gv.;(x)*(x)«(^)*(^)/(« j )[i-/(«)] ( (l ,_ (et _ £ 1 . )/g))+i0+ - („ + ( a y )M ) +ffl + ) 

(S17) 

In Fig. [5] we show the oscillation spectrum at finite temperature in the TG case obtained with the real-time 
evolution presented in Section IV.A above and with the linear response theory approach. The results obtained with 
both methods agree quite well, which is all the more remarkable given the finite simulation time used for the real-time 
evolution calculation. 


|Ke[dxc M («j)]| 



u/u h 


FIG. 5. (Color online). Fourier transform of the center of mass oscillation |5fte[<kc C M(<^)] | as a function of frequency, in the 
TG limit obtained with the Fourier transform of the real-time evolution (black solid line) and from linear response theory (red 
vertical lines with dots). 


Density profiles 


In Fig. [6] we show the density 


profiles in the presence and absence of the barrier for different interaction regimes. 



FIG. 6. Ground state density for different interparticle interaction regimes and Uo/Eh^h = 0,1 (dashed and solid line respec¬ 
tively). In the TG case the barrier creates a notch in the density, in correspondence of a minimum or maximum, depending on 
whether the number of particles is even or odd. The grey line is a sketch of the potential in dimensionless units and divided 
by a factor of 50. 
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Scaling of the dipole frequency shift with number of particles 


In Fig. [7] we show the scaling behavior of the dipole frequency as a function of the number of particles TV, as 
obtained from the exact calculation in the TG regime. For sufficiently large numbers of particles we find a perfect 
agreement with the |cjd — co>h| °c l/\/TV power law scaling predicted with the LL technique. 



N 


FIG. 7. (Color online). Scaling of ua as a function of the number of particles N in the TG regime. Different colors correspond 
to A = Uo/Ehlh = 0.2,1, and 5 (orange, red, and brown points respectively). For N sufficiently large (N > 10) the scaling 
follows the power law |<x?d — cdh| oc 1 /\/N (black solid lines). 
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